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It is well known that anyons in topologically ordered 
phases can carry symmetry quantum numbers that are 
quantized to fractional values. In the celebrated exam¬ 
ple of fractional quantum Hall states, Laughlin quasi¬ 
particles carry fractional charge—the quantum number 
of the U( 1) symmetry jTj. In recent years, great progress 
has been made in understanding the interplay between 
symmetry and fractionalization in other topologically or¬ 
dered states. In particular, topological spin liquids ex¬ 
hibit a more subtle kind of symmetry fractionalization, 
associated with the crystal symmetry of the underlying 
lattice instead of internal symmetries W- While some 
aspects of it have been studied for quite a while, crys¬ 
tal symmetry fractionalization has now received renewed 
attention, due to an increased interest in the role of crys¬ 
tal symmetry in topological phases of matter. This topic 
is also becoming timely in view of strong numerical evi¬ 
dence for spin liquids on kagome lattice found in the last 
few years [Ml- In order to fully pin down the topolog¬ 
ical nature of the numerically found spin liquid liquid, 
the complete pattern of crystal symmetry fractionaliza¬ 
tion needs to be determined. 

In this work, we offer a new perspective on crystal 
symmetry fractionalization in Z 2 spin liquids. We find 
that the nontrivial way that crystal symmetry acts on 
an individual anyon is directly related to the symmetry 
representation of the topologically ordered ground states, 
as labeled by the crystal momentum and parity of many- 
body wave functions. Given that states with different 
symmetry labels cannot be adiabatically connected, our 
finding immediately makes it clear that the classification 
of spin liquids is refined and enriched by taking into ac¬ 
count crystal symmetries HD]. Our theoretical result also 
provides a straightforward method to classify and detect 
different spin liquids in numerical studies. As a concrete 
example, we demonstrate that our method can be used to 


easily distinguish various Z 2 spin liquids on the kagome 

lattice HIKED]. 

We begin by briefly reviewing what is known about 
crystal symmetry fractionalization in Z 2 spin liquids, and 
setting up the terminology for our work. A Z 2 spin liq¬ 
uid H2 ITS] supports three types of anyon excitations: 
bosonic spinons, fermionic spinons, and visons. As a 
defining property of topological excitations, anyons of 
each type can only be created in pairs. This property 
makes symmetry fractionalization possible. This can be 
understood by considering a many-body excited state 
containing two identical anyons that are spatially sep¬ 
arated 13 HD]. Intuitively speaking, the action of sym¬ 
metry on this excited state can then be factorized into 
a product of two independent symmetry actions on the 
anyons. While the action on a physical state is necessar¬ 
ily described by a linear representation of the symmetry 
group denoted by G, the action on a single anyon is now 
allowed to form a projective representation G, such that 
the tensor product G ® G is a linear representation of 
g HMD]. The projective representation G, which has a 
different group algebra than G, can be regarded as the 
“square root” of G. In this sense, symmetry action on 
anyons can be called “fractionalized.” Throughout this 
Rapid Communication, a tilde symbol (~) placed over 
a symmetry operation means that it acts on an anyon; 
otherwise it acts on a physical wave function. 

To sharpen the intuitive argument stated above and 
give a precise definition of symmetry fractionalization is 
a nontrivial task that requires great care. The main diffi¬ 
culty is that symmetry operations should in principle be 
performed on a single anyon, and yet any physical wave 
function necessarily contains an even number of them. 
(We note that crystal symmetry fractionalization may 
have implications for excitation spectra that can be de¬ 
tected To overcome this difficulty, we take a 
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different approach and give a precise and operational def¬ 
inition of crystal symmetry fractionalization by relating 
it to the linear symmetry representations of many-body 
wave functions. 

Crystal symmetries of a given lattice form a space 
group G generated by translation, rotation and reflection. 
Any allowed projective representation of G, denoted by 
G, can be specified by its modified group algebra as com¬ 
pared to G[ig. First, two commuting operations A and 
Y in G, XY = YX, can become anti-commuting in G, 
XY = — YX. This will be referred to as commutation 
relation fractionalization. Second, an identity of X n = 1 
in G can become X n = — 1 in G. This will be referred to 
as quantum number fractionalization. 

As the main result of this work, we find that the 
symmetry representation of ground states is a diagno¬ 
sis of crystal symmetry fractionalization in Z 2 spin liq¬ 
uids. First, the commutation relation between a trans¬ 
lation operation Tj and another symmetry operation A, 
TiX = ±X7j, can be determined from the difference 
between eigenvalues of X for ground states in different 
topological sectors on a torus geometry with an odd num¬ 
ber of unit cells in the direction of Tj. Second, for an 
order-two symmetry operation X 2 = 1, X 2 = ±1 can 
be determined from the parity eigenvalue of X acting on 
ground states on a torus with 4n + 2 sites. We find it 
remarkable that the fractionalized symmetry property of 
anyons is simply and directly encoded in the symmetry 
of ground states. Some technical details of our derivation 
are available in the Supplemental Material 1 125] , 

Fractionalized commutation relation. In the presence 
of a fractionalization in the commutation relation be¬ 
tween T\ and another symmetry operation A', ground 
states in different topological sectors have different eigen¬ 
values of A on a torus with odd number of unit cells in 
the Ti direction [5]. To extract the fractionalization of 
different anyons, we find it crucial to choose a basis ac¬ 
cording to the anyon flux along the direction of T\. 

On a torus, a gapped Z 2 spin liquid has a fourfold 
ground state degeneracy, which is protected by its in¬ 
trinsic topological order. A basis of these four ground 
states can be chosen such that each state carries a differ¬ 
ent anyon flux going through the torus in the direction of 
Ti , which can be diagnosed by a Wilson loop operator in 
the direction of T 2 |T5] . There are four types of such flux; 
each corresponds to one type of anyon excitations in the 
toric code topological order. Therefore the four ground 
states can be labeled as |G a ), where a denotes the type 
of anyon. In this paper we denote the four types of anyon 
excitations in the Z 2 spin liquid state as a = 1, 6, /, and 
v, standing for the trivial particle, the bosonic spinon, 
the fermionic spinon, and the vison, respectively. 

Using this basis, the ratio between parity eigenvalues 
of two ground states can be calculated by considering the 
Berry phase picked up through the following operations 


acting on |Gi), 

A- 1 (/°)- 1 Ar|G 1 ) = e iA$ |G 1 ) ! (1) 

where f a denotes the operation of moving one a anyon 
across the torus in the direction of T\ and it maps |Gi) 
to / a |Gi) = |G a ) [24]. This Berry phase A$ can be 
obtained from the ground state A-symmetry representa¬ 
tions as the following, 

\Gi) A |G„) ^ \ a x \G a ) 

A ^| Gl) 2L1+ Al-(Ai-)- 1 |G 1 ), 

where A^ denotes the parity eigenvalue of A acting on 
|G„>: A|G„) =\ a x \G a ). 

On the other hand, the same Berry phase can be ob¬ 
tained using the projective crystal symmetry represen¬ 
tation of the a anyon. Starting from the ground state 
|Gi), the operation f a creates two anyons locally and 
moves one across the torus along Tj, which is equivalent 
to acting T\ on one of the two anyons n\ times, so the 
end state can be expressed as a®T? 1 a, where ® denotes 
anyon fusion and n\ is the number of unit cells in the 
direction of Ti , which is an odd number according to our 
setup. Then A acts on both anyons and maps the state 
to Xa® XT ™ 1 a. The rest of the actions can be calculated 
similarly, 

1 = a®a ^ a® fj ni a ^ A a ® XT? 1 a 

{fa) ---> Xa ® T? ni XT? 1 a a ® X~ 1 f? ni XT ? 1 a. 

Hence after the series of operations one anyon is changed 
into X~ 1 T? ni XT? 1 a = (r£) ni a, where r ? = ±1 de¬ 
notes the fractionalization of the commutation relation, 
TjAa = t?X1\ a. This can be simplified as t x because 
ni is odd and (r x ) 2 = 1. Comparing this result with 
Eq. ([ 2 ]) we obtain the following relation between commu¬ 
tation relation fractionalization and ground state parity 
eigenvalues, 

T x = ^x/^x- ( 3 ) 

Fractionalized quantum number. The action of an 
order-two crystal symmetry A on an anyon is fractional¬ 
ized if acting A twice on a single anyon yields A 2 = — 1. 
To detect such fractionalized quantum number, we act 
A once on an excited state containing two anyons, whose 
positions are swapped by A [25]. Specifically, A maps an 
anyon at a site i to another anyon at the image site A (i) 
and vice versa. 

The symmetry action on an anyon is accompanied by 
additional gauge transformations [2], 

Aoj = Uia X (i), Xa X (i) = U x ^ai. (4) 

Therefore, acting A twice on the anyon a, leaves anyon at 
its original position but yields a factor A 2 a = U x ^Uia. 
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Alternatively, if one perform the operation X once on a 
physical wave function |’F) that contains a pair of anyons 
at i and X(i), the same phase factors Ui and Uxu) are 
collected from each anyon, so that the pair of spinons 
acquires the same total phase of U X (i)Ui- Importantly, 
if the anyon under consideration is a fermion, there is 
an additional statistical sign due to the exchange of two 
fermions under X. To summarize, when a bosonic anyon 
carries a fractionalized quantum number X 2 = — 1, the 
parity eigenvalue of an excited state I'll) containing a pair 
of such anyons is opposite to that of the ground state, 
from which |\H) is created. For fermionic anyons, X 2 = 
— 1 implies that the parity eigenvalues of |\H) and |G) are 
identical. 

We now show that detecting X 2 = ±1 for spinons can 
be further simplified when the spin liquid is constructed 
from parton methods (using either the Schwinger-boson 
or Abrikosov-fermion approach), which put exactly one 
spinon on every lattice site. Specifically, we choose a 
lattice geometry with 4ra + 2 sites, so that the ground 
state contains an odd number of pairs of spinons. The 
parity eigenvalue of such a ground state is then equal 
to the parity of a single pair of spinons, which in turn 
detects X 2 = — 1 as described above. Importantly, the 
contribution to the parity eigenvalue from the fractional 
quantum number is independent of the topological sec¬ 
tor of the ground state, in contract to the previous case 
involving the commutation relation between A' and T\. 
This will be demonstrated with concrete examples below. 

Z 2 spin liquids on kagome lattices. We now apply our 
method of detecting crystal symmetry fractionalization 
to Z 2 spin liquid states on the kagome lattice. Recent 
numerical studies of the spin-1 Heisenberg model on the 
kagome lattice |5HS] have found strong evidence for a 
gapped spin liquid state, likely with a Z 2 topological or¬ 
der (see however, Ref. (21j). On the other hand, various 
types of Z 2 spin liquids on the kagome lattice that dif¬ 
fer in symmetry properties have been theoretically con¬ 
structed using parton methods in early studies inHsi- 
In what follows, we will connect the numerical findings 
to theoretical constructions and show how to determine 
which of the Z 2 spin liquid states theoretically proposed 
so far is consistent with the ground state of the Heisen¬ 
berg model on the kagome lattice. 

To start, we quickly describe the parton construction of 
various Z 2 spin liquid states, paying particular attention 
to the role of crystal symmetry. Parton constructions 
postulate that the low-energy dynamics of the spin liq¬ 
uid phase is described by gapped spinons (which carry Z 2 
gauge charge) interacting with Z 2 gauge fields. Depend¬ 
ing on the Z 2 background flux patterns, the spinon ex¬ 
hibits different crystal symmetry fractionalizations, thus 
leading to distinct spin liquid states. This is because in 
the presence of gauge flux, crystal symmetries acting on a 
spinon involve additional Z 2 gauge transformations. For 
example, when there is a n flux within a unit cell, trans¬ 


lations of spinons correspond to the magnetic translation 
group with the property T 1 T 2 = —T 2 T 1 , resulting in a 
fractionalized commutation relation. 

In the parton construction, spin liquids with differ¬ 
ent Z 2 background flux patterns are classified using the 
projective symmetry group (PSG) analysis invented by 
Wen [2], from which we can derive crystal symmetry 
fractionalization of bosonic and fermionic spinons. Be¬ 
low we derive the crystal symmetry fractionalization for 
Z 2 spin liquids on the kagome lattice that were con¬ 
structed using the Schwinger-boson approach in previous 
works [111 fl2i . The PSG analysis by Wang and Vish- 
wanath [12] has found that there are four spin liquid 
states with different Z 2 flux patterns, which are adiabat- 
ically connected to nearest-neighbor resonating-valence- 
bond states and therefore have better variational energy 
than other states. Hence we use them as examples to 
demonstrate our method of detecting crystal symmetry 
fractionalization. These four states are labeled by three 
Z 2 variables {pi,P 2 ,P 3 )', for the sake of completeness, this 
terminology from Ref. [l2]is reviewed in Sec. I of the Sup¬ 
plemental Material. 



FIG. 1. The definition of crystal symmetry operations. 7i,2, 
o, and R * label the symmetry operations of translation, mir¬ 
ror reflection, and sixfold rotation, respectively. 

The kagome lattice has three independent symmetry 
operations: Tj, 2 , cr, and R n / 3 , which denote translation, 
mirror reflection and rotation, respectively, and their def¬ 
inition is shown in Fig. [l] A straightforward translation 
of terminology shows that (pi,P 2 ,P 3 ) in the PSG analysis 
directly yields the crystal symmetry fractionalizations of 
the bosonic spinon excitations (denoted by b), listed in 
the first row of Table [Tj Here, ry 2 = ±1 labels the frac¬ 
tionalization of the commutation relation between T) and 
T 2 , defined by TjTj = Ty, TjXj. Likewise, T a is defined by 
Xj<r = r a T\(7, and the fractionalization of commutation 
relation between Tj and the twofold rotation R n = R^/ 3 

takes the form of TjfG = tr^R^T^ . All three r’s are 
equal for the four spin liquid states we consider. In ad¬ 
dition, quantum number fractionalizations d 2 = ±1 and 
R 6 = ±1 are listed in the last two columns of Table II. 

TV j o 

A limitation of the previous PSG analysis is that 
it is tied to the Schwinger-boson formalism and hence 
only gives the crystal symmetry fractionalization of the 
bosonic spinons. The vison excitation in all four states 
has the same crystal symmetry fractionalization HUH: 
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TABLE I. Crystal symmetry fractionalizations of different 
anyon excitations, b, /, and v denotes the bosonic spinon, 
the fermionic spinon and the vison, respectively. 


Anyon 

TTi = 

a 2 

06 p 2 

ix 2L fl'TT 

b 

(- 1) P1 

(_ 1 )P 2 

(_1)P1+P3 

f 

(-1) P1+1 

(_1)P2 + 1 

(_1)P1+P3 + 1 

V 

-1 

+ 1 

+1 


t t 2 = T a = tr^ = — 1 and a 2 = R 2 = 1. This result 
can be simply obtained from the charge-flux duality: in 
a spin liquid state with odd number of spinons per unit 
cell as is the case for the kagome lattice, the vison always 
sees a n flux per unit cell because a spinon is a n flux to 
a vison. As a result, the vison always has the property 
TiT 2 = -f 2 f,. 


We now use the method of flux attachment to derive 
the crystal symmetry fractionalization for the fermonic 
spinon, which is equivalent to a composite of a bosonic 
spinon and a vision—the latter is a n flux to the former. 
Compared to a bosonic spinon, the fermionic spinon al¬ 
ways sees an extra n flux per unit cell due to the vison 
attached to it; hence Tj- 2 = — . Similarly, the differ¬ 

ence in R 2 between bosonic and fermionic spinons follows 
from the fact that the attachment of a n flux changes the 
angular momentum between integer and half-integer val¬ 
ues [25j. These results can also be derived using general 
methods described in Sec I of the Supplemental Material, 
and are summarized in the second row of Table [I] [29]. 


TABLE II. Correspondence between Schwinger-boson and 
Abrikosov-fermion constructions, pi labels the PSG solutions 
of Schwinger-boson construction [T2]. The Q\ = ±Q 2 labels 
are used by Sachdev ill], and the labels in the last column 
are used by Lu et al. US¬ 


(Pl,P2,P3) 

Label in Ref. [IT] 

Label in Ref. OGO 

(0,0,1) 

Q 1 = —Q 2 

Z 2 [0,7r]a 

(0,1,0) 

Qi = Q 2 

Z 2 [0, 7t]/3 

( 1 , 0 , 1 ) 


Z 2 [0,0]B 

( 1 , 1 , 0 ) 


Z 2 [0,0]A 


As an independent check of the above results, we find 
that the above four spin liquid states constructed from 
the Schwinger-boson approach can be equivalently de¬ 
scribed using the Abrikosov-fermion approach. To es¬ 
tablish the mapping between the two parton construc¬ 
tions, we match the ground states of spin liquids in the 
nearest-neighbor resonating-valence-bond limit, given by 
the Gutzwiller projection on the corresponding parton 


wavefunctions in the two constructions: 


= P G exp 


= P G exp 


£,ij e <xpb\ a b*. 




(ij) 


Y. Cij e apfi a fjp 


(«> 


| 0 >, 

| 0 >. 


(5) 

( 6 ) 


Here Gj and Qj are antisymmetric and symmetric scalars 
on the nearest-neighbor bonds (ij) used in Schwinger- 
boson and Abrikosov-fermion constructions, respectively. 
The values of Gj or Qj on bonds of the kagome lat¬ 
tice are given by the corresponding PSG analysis. By 
equating Eq. <[5j) and §• we find explicitly the map¬ 
ping between and Qj for each of the four spin liq¬ 
uid states [301 EH- As a byproduct, this mapping es¬ 
tablishes the correspondence between the notation of the 
Schwinger-boson mm formalism (pi,P 2 ,P 3 ) with that 
of the Abrikosov-fermion }13l formalism, as shown in Ta¬ 
ble [IT] From this mapping, we have confirmed the results 
in Table [T] (see Sec II of the Supplementary Material for 
details). 


TABLE III. Ratios between A-symmetry parity eigenvalues 
of | G a ) and |Gi}. The results depend on the commutation 
relation fractionalization r|> 2 = — t£ 2 = (—1) P1 , and the ratios 
are the same for A = T 2 , t and R n . 


Pi 

\ b x/Xx 

Xk/Xx 


0 

+1 

-1 

-1 

1 

-1 

-1 

+1 


Having derived the crystal symmetry fractionalization 
of anyons for all four spin liquids on the kagome lattice, 
we can now determine the symmetry representations of 
the ground states for each spin liquid, by using the gen¬ 
eral relation between the two as described in the first 
part of this work. First, based on Eq. ([ 3 ]) and the results 
of commutation relation fractionalizations summarized in 
Table [T] the ratios between symmetry eigenvalues of dif¬ 
ferent ground states are determined and the results are 
summarized in Table m 

Second, if all other aspects of the symmetry fraction¬ 
alization are the same, a difference in the quantum num¬ 
ber fractionalization results in a uniform parity change 
in the symmetry representation of ground states in all 
topological sectors. This relation can be obtained by ex¬ 
plicitly calculating the parity eigenvalues of the model 
wave functions in Eq. ((5) and Eq. (|6|. The results are 
summarized in Table Mil and the details of the deriva¬ 
tion is given in Sec. Ill of the Supplemental Material. 
These results are determined from projected mean field 
wave functions but they also apply to general wave func¬ 
tions in the same topologically ordered phase, because 
the crystal symmetry representations are invariant when 
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the state is smoothly deformed without breaking crystal 
symmetries. 


TABLE IV. Crystal symmetry representations of ground 
states in different topological sectors on a torus with odd- 
by-(4n + 2 ) unit cells. \G a ) denotes the ground state with an 
anyon flux a in the direction of T\, where a = 1, b, v, and / 
denotes the trivial anyon, bosonic spinon, the vison and the 
fermionic spinon, respectively. 


X 

\Gi) 

|G fc ) 

1 Gv) 

1 G f ) 

t 2 

1 

(-IP 

-1 

(_l)Pl + i 

a 

(-1) P2 

(_ i)Pi+P 2 

(_l)P2 + i 

(_1)P1+P2+! 

Rw 

(i r 

(_ l)Pi+P3 

(_1)P3 + ! 

(1)P1+P3+ 1 


Summarizing the above results, we see that the PSG 
parameters p 2 and p 3 determines the parity eigenvalues of 
the ground state |Gi), and then using pi the parity eigen¬ 
values of the other sectors are also determined. Hence we 
can obtain from (pi,P 2 ,P 3 ) the crystal symmetry repre¬ 
sentations of all topological sectors on a torus with odd- 


by-(4n + 2) unit cells, as summarized in Table IV 


Conclusions. In this work we propose a method to 
detect crystal symmetry fractionalizations from the crys¬ 
tal symmetry representations of the ground states. On 
a torus with 4n + 2 sites, the ratio between symmetry 
parity eigenvalues of ground states in different topologi¬ 
cal sectors detects the fractionalization of commutation 
relation with a translation symmetry, and a uniform sign 
change in all sectors detects the quantum number frac¬ 
tionalization. 

Our method can be applied to study the nature of the 
topological order of the Z 2 spin liquid states obtained 
in numerical studies. Particularly using the infinite- 
size density matrix renormalization group (DMRG) 
method [32] or an infinite projected entangled pair states 
(PEPS) ansatz [33] the ground states in different topo¬ 
logical orders can be obtained on an infinite cylinder and 
labeled with the anyon type by calculating the modu¬ 
lar matrices [H US Eg. To study both commutation 
relation fractionalization and quantum number fraction¬ 
alization, we suggest using a cylinder (4 n + 2) unit cells 
wide. Then the crystal symmetry fractionalization stud¬ 
ied in this work can be determined by putting the system 
on a torus with a length of an odd number of unit cells 
and examining the crystal symmetry representation of 
topologically degenerate ground states. For the Z 2 spin 
liquid states discussed above, the results are shown in 
Table HVl 

In this work we explicitly derived the crystal symme¬ 
try quantum number of ground states and the fractional 
symmetry quantum number of anyons for four Z2 spin 
liquid states on the kagome lattice m ■ Our method can 
be used straightforwardly to study additional spin liquid 
states that have been theoretically proposed [u uni [ 13 ]. 
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Supplemental Material 


In this supplemental material we provide some technical details of the derivations outlined 
in the main text. 


I. CRYSTAL SYMMETRY FRACTIONALIZATIONS OF BOSONIC AND 

FERMIONIC SPINONS. 


In this section we derive the relation between crystal symmetry fractionalizations of 
bosonic and fermionic spinons. In a Z 2 spin liquid states, both of the three types of topolog¬ 
ical excitations (bosonic spinon, fermionic spinon and vison) carry projective representations 
of the crystal symmetry group, and different anyon carries a projective representation with 
different crystal symmetry fractionalizations. However, the crystal symmetry fractionaliza¬ 
tions of different anyons are not independent, and particularly the symmetry fractionaliza- 
tion of fermionic spinon can be derived from the ones of the bosonic spinon and vison, as the 
former is a bound state of the latter [1], For Z 2 spin liquids with odd number of spinons per 
unit cell, the symmetry fractionalization of the vison is fixed [2, 3]. Therefore the symmetry 
fractionalization of the fermionic spinon can be uniquely determined from the one of the 
bosonic spinon. In this section we focus on the following three crystal symmetry fractional¬ 
izations on the kagorne lattice: the commutation relation fractionalization between Tj and 
X 2 , and the quantum number fractionalization of a 2 and R 2 . These results are used in the 
main text to obtain Table I. 

Firstly, the fractionalization of commutation relation between Tj and X 2 of fermion, r^, 
can be derived from the corresponding fractionalizations of the bosonic spinon and of the 
vison [1], 


t t 2 = t t 2 t t 2 - 


( 1 ) 


For Z 2 spin liquids with odd number of spinons per unit cell, t^ 0 = —1 [2], Hence we always 
have Ty 2 = —Tj - 2 [3]. As explained in the main text this can be understood intuitively using 
flux attachment. 

Next, we consider the reflection symmetry a. We find as a bound state of a bosonic 
spinon and a vison, the fractionalized quantum number of a fermion is the product of the 
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quantum numbers of a bosonic spinon and that of a vison, times an additional twist factor 
t = -1 

( 2 ) 

where a a denotes the projective representation of a acting on anyon type a. The twist factor 
t — — 1 can be calculated by examining the reflection parity eigenvalue of a wave function |T) 
containing two fermionic spinon excitations located at two sites related by cr. As discussed 
in the main text, if aj = +1, the reflection parity eigenvalues of |T) and of the ground state 
\G) are opposite because of the Fermi statistical sign of exchanging two fermions, and if 
aj = — 1 the two eigenvalues are the same. On the other hand, IT) can also be viewed as a 
wave function containing two bosonic spinons and two visons. Therefore the ratio between 
reflection parity eigenvalues of |T) and of |G) is also given by dfcd;. Comparing these two 
results we conclude that the twist factor t = — 1 in Eq. (2). For the Z 2 spin liquid states we 
consider, the vison always has of = +1 [2], Hence the fermionic spinon and bosonic spinon 
always have opposite reflection quantum number fractionalization aj = — of. 

Our conclusion of t = —1 should be contrasted with the incorrect result of the twist factor 
obtained in Ref. 1. We note that the phase of the two-fermion wave function depends on the 
ordering of the two fermions, and after mirror reflection the exchange of the two fermions 
gives the Fermi statistical sign of —1. This minus sign is the key in the above derivation of t = 
— 1. As a independent check, in Sec. 11 we explicitly derive the mapping between Schwinger- 
boson and Abrikosov-fermion constructions, which gives the same mapping between the 
crystal symmetry fractionalizations of the bosonic and fermionic spinons. 

We also note that a different operational definition was used by Lu et al. [3], which lead 
to the identical result. In their definition they considered two anyons, whose positions are 
unchanged under mirror reflection, in contrast to our definition where the position of the 
anyons are exchanged by mirror reflection. Their definition is more involved to study when 
the translation operator (which connects the two positions of the anyons) does not commute 
with the reflection. 

Similarly for the inversion symmetry the twist factor t is also —1 because of the Fermi 
statistical sign, so Rjj = —Rj b Rj v . Again for the Z 2 spin liquids on the kagome lattice we 
always have Rj v = +1. This twist factor t — — 1 was also obtained using different methods 
by Essin and Hermele [1]. Therefore the fermionic spinon and bosonic spinon have opposite 
inversion quantum number fractionalization R‘j j = —Rj b . 
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II. ABRIKOSOV-FERMION REPRESENTATION OF SCHWINGER-BOSON 

SPIN LIQUID STATES 


In this section we discuss the Abrikosov-fermion representation of the nearest-neighbor 
Schwinger-boson spin liquid states. In general these two ways of constructing spin liquid 
states can describe several common Z 2 spin liquid states [3], and particularly spin liquid 
states with pairings on nearest-neighbor bonds only can be described by both of the two 
constructions [4, 5]. 

We begin with a brief review of the four bosonic PSG solutions discussed in the main 
text, which were studied in details by Wang and Vishwanath [6]. These four PSG solutions 
are labeled by three Z 2 variables (pi,p 2 ,p:i), which parameterizes the projective crystal 
symmetry representation carried by the bosonic spinon operators, 


<hi ( r ) = 0 , 


(3a) 


4>T 2 {r) =Pi7rr 1 , 


(3b) 




71 

= P2~ + PlKT 2 , 

vr r 2 (r 2 - 1 + 2rP 

= P*2 + P - 2 - 


(3c) 

(3d) 


where (j>x{?) denotes the additional U(l) phase acquired by the bosonic operator after the 
symmetry operation X. In other words, the projective symmetry representation of X carried 
by b ia is 

X : b ia -A e^xMbxva. (4) 

The PSG solution in Eq. (3) determines both the crystal symmetry fractionalization 
of the bosonic spinon and the pattern of in the mean held wave function in Eq. (6) 
of the main text. On one hand, the PSG solution describes the projective representation 
of the crystal symmetry group carried by the bosonic spinons, and from it their crystal 
symmetry fractionalization listed in Table I can be derived. On the other hand, it also 
describes the pattern of ^ such that the wave function in Eq. (6) is invariant under symmetry 
transformations. The patterns of for the four different PSG solutions are plotted in Fig. 1, 
where on each bond the value of ^ l3 = ±1 is represented by an arrow because ^ = ~^ 3 i- 
These patterns are the same as the patterns of the nearest-neighbour pairing terms listed in 
Ref. 6. 
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FIG. 1. Pairing amplitude of Schwinger-boson mean field constructions with four different 
PSG solutions. The arrows represent the value of Gj = ±1, where an arrow from i to j indicates 
£ij = ~~£ji = Tl. 

TABLE I. Summary of PSG solutions for the Schwinger-boson and Abrikosov-fermion construction 
of the four Z 2 spin liquid state. 


(p b l,P2,Pl) 

{p{,pLpI) 

(1,1,0) 

(0,0,0) 

(0,1,0) 

(1,0,0) 

(1,0,1) 

(0,1,1) 

(0,0,1) 

(1,1,1) 


Starting from a specific Schwinger-boson PSG solution, we derive the PSG solution of 
the corresponding Abrikosov-fermion PSG by first working out the pattern of Q 3 using 
the correspondence between ^ and Qj given by Ref. 5, and then identify the compatible 






FIG. 2. Pairing amplitude Qj of Abrikosov-fermion mean field constructions with four different 
PSG solutions. A thin (thick) bond represents Qj = +1 (-1), respectively. 
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PSG solutions, which are described by the other four combinations of (pi,P 2 ,P 3 ) that does 
not describe a nearest-neighbor Schwinger-boson ansatz. These Abrikosov-fermion PSG 
solutions are a subset of the ones given by the general classification in Ref. 7. Here we 
denote the PSG parameters of the Schwinger-boson and Abrikosov-fermion ansatz by p\ and 
p{, respectively. For each Schwinger-boson state labeled by (p\,p\,pX)i the corresponding 
Abrikosov-fermion PSG solution (p {, pf 2 , p^) and the label as used in Ref. 7 is listed in Table I. 

In these four Abrikosov-fermion PSG solutions, the fermionic spinons transform under 
crystal symmetry operations with an additional U(l) phase described by Eq. (3), and these 
solutions are a subset of the possible PSG solutions, with the general SU(2) gauge transfor¬ 
mations restricted to the U(l) subgroup. Hence these four states are included in the general 
classification given by Lu et al. [7]. We can label these states with the notation used in Ref. 7 
by calculating the topological invariants of symmetry fractionalizations listed in Ref. 3. 


TABLE II. Summary of PSG solutions and exactly solvable Hamiltonians of the four different Z 2 
state realizable in the quantum dimer model. 


Symmetry fractionalizations 

(Pi,P2>pI) 

SU(2) PSG 

T-'T- l T 2 Ti 

(-1 )p( 

P 12 

a 2 

(-1)^ 

PaPcrCe 


(_1)P2+P3 

Per 

(^f) 6 

(_1)P{+P3 

PC 6 


("I ) Pi 

PaTPC 6 T 

S3 

S3 

GO |-1 

(-1)^3 

PC 6 T 


TABLE III. Summary of PSG solutions and exactly solvable Hamiltonians of the four different Z 2 
state realizable in the quantum dimer model. 


(Pl:P2iPi) 

P12 

Pa- 

PaT 

PaCe 

PC 6 T 

PCs 

SU(2) PSG Label 

(0,0,0) 

+1 

+1 

+ 1 

+ 1 

+ 1 

+ 1 

Z 2 [0,0]A 

(1,0,0) 

-1 

+1 

+ 1 

+ 1 

+ 1 

-1 

Z 2 [0, tt\/3 

(0,1,1) 

+1 

+1 

+ 1 

-1 

-1 

-1 

z 2 [0,0]b 

(1,1,1) 

-1 

+1 

+ 1 

-1 

-1 

+ 1 

Z 2 [0,7r]a 


Using the PSG solutions in Eq. (3), six symmetry fractionalization invariants are ealeu- 
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lated and listed in Table II. These results are then compared to the same invariants expressed 
using the Z 2 labels t)a for generic SU(2) PSG solutions[3, 7]. Using these relations we can 
cauculate the labels tja for the four PSG solutions labeled by p{ and the results are listed 
in Table III. Finally these states can be labeled by the notations defined by Lu et al. [7] 
by comparing the results of r )a to Table III in Ref. 3, and these labels are listed in the last 
column of Table III. 

In summary, starting from a nearest-neighbor pairing Schwinger-boson state with a cer¬ 
tain PSG solution, we can first derive the corresponding fermionic PSG solution labeled by 
{p{iP{iP{) using the mapping in Table I, then label the PSG solution using the notation in 
Ref 7 using the mapping in Table III. These results are presented in Table II in the main 
text. 


III. PARITY EIGENVALUES OF THE GROUND STATE WAVE FUNCTIONS. 

In this section we calculate the a and R n parity eigenvalues of the ground state wave 
functions, which is used to detect the fractionalized quantum numbers of a 2 and R as 
explained in the main text. We first obtain the results using the model wave functions in 
Eq. (6) and Eq. (7) in the main text, then argue that this result also applies to generic wave 
functions that belong to the same SET phase. 

Here we put our system on a torus with 4n + 2 sites, so there must be an odd number of 
unit cells along one direction, and as shown in Table III in the main text different ground 
states have different parity eigenvalues of a and R n . Hence we also need to determine the 
topological sector the mean field wave functions belong to. Using the sign structure in 
the mean field wave functions one can show that [8] the Schwinger-boson wave function in 
Eq. (6) in the main text is a superposition of |Gi) and |Gft), while the Abrikosov-fermion 
wave function in Eq. (7) in the main text is a superposition of (S',,) and |G^). Using this 
result it is easy to identify the ground state crystal symmetry representations from the 
projected mean field wave functions. 

For the two states with p\ = 0, the crystal symmetry representations can be determined 
using the Schwinger-boson wave function. In this case their PSG solutions in Eq. (3) are 
simplified to that after the action of o or R n the spinon operator bi a acquires a phase of 
p\ | and respectively [6]. Therefore the wave function, containing 4n + 2 spinons after 
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the Gutzwiller projection, acquires a phase of p 2 ir and pgTr, respectively. Since this wave 
function is a superposition of G'i) and Gb), we conclude that their parity eigenvalues under 
a and R^ is (—1) P2 and (—l) p y respectively. 

For the two states with p\ — 1, we need to use the Abrikosov-fermion wave function. (In 
this case the Schwinger-boson wave function is not symmetric because it is a superposition of 
|Gi) and \Gb), which has different crystal symmetry representations in this certain geometry 
according to Table III in the main text.) Similarly in these two states the PSG solutions also 
simplies to that after the action of cr or R n the spinon operator f ia acquires a phse of p 2 § and 
pi |, respectively. However in this case in addition to this phase from each spinon, there is an 
additional sign from the reordering of fermion operators in the Gutzwiller projection operator 
after the symmetry transformation, which is worked out in Sec. IV. Combing these two phase 
factors we conclude that the parity eigenvalues of the Abrikosov-fermion wave function, 
which is a superposition of \G V ) and | Gb), is (—1) P2 = (—1)P 2+1 and (—1)P3+ 1 = (—ijp-H 1 
under a and R n , respectively (we used the mapping between p\ and p{ in Table I). Hence 
according to the results in Table III in the main text the parity eigenvalues of |Gi) are 
(-l ) p 2 and (-l) p 3. 


IV. FERMION REORDERING SIGNS. 

In this section we compute the fermion operator reordering signs that appears in the 
definition of the Gutzwiller projection operator for the Abrikosov-fermion mean field state 
in Eq. (7) in the main text. The projection operator Pq contain a product of fermion 
operators on each lattice sites, and after a crystal symmetry transformation the ordering of 
the sites is changed, so the aforementioned wave function can acquire an additional minus 
sign due to the reordering of fermion operators. This sign qr x , where X denotes the crystal 
symmetry operation, is computed here, for a system with 4n + 2 sites. 

For the translation symmetry operation Ti, we can order the fermion operators first from 
left to right in each row in the direction of 7\, then from top to bottom for each row in the 
direction of T 2 . After translation, in each row we need to move one fermion operator from 
the beginning to the end. Assuming there are n i and n 2 unit cells in the direction of Tj and 
T 2 , respectively, the reordering of the fermion operators gives a sign of (—l) ni+1 for each 
row, and a total sign of (—l) n2 ( ni + 1 ) — (—i) n2 because the total number of unit cells n\ n 2 



is odd. Since we have an odd number of sites per unit cell and each sublattice is giving the 
same contribution, the total phase is ip* T = t^tt. Similarly we have (pj, = ni7r. Note that 
since ri\ri 2 = 4 n + 2, one of these two factors is n and the other is zero. 

For the inversion symmetry operation R w , we can order the fermion operators in a way 
that the two sites related by R n appear together: 

w)p = {flA) (■ ftA ) • • • (4„ +3 4 +3 ) i°)’ ( 1 2 3 4 5 6 7 8 ) 

where // = fLi an( l q labels the lattice site that is related by R n to site i. Since there 
are in total 4n + 2 unit cells in the system, there are in total 6 n + 3 pairs of sites. After the 
operation of R n , we need to exchange the two fermion operators within each pair to restore 
the fermion order, and each pair contributes a minus sign. Since the total number of pairs 
is odd, we conclude tp] f Rr = n. 

For the mirror reflection symmetry cr, we also order the fermions such that two sites 
related by cr appear together. However, as shown in Fig. 1 in the main text, there are also 
sites that are cut by the mirror plane and stays invariant under mirror reflection, and we 
put them in the beginning of the product, as the following, 

Wv -4-4. (44) (44) ■ ■' (44) i°>> («) 

where j\, ■ ■. j p denotes the sites that are on the mirror axis cr. The total number of such 
sites p equals to the number of unit cells along the direction of cr. In our setup, the number 
of sites on the mirror plane cr always equal to two modular four. Hence the number of pairs q 
in Eq. (6) must be an even number. Similar to the previous discussion we conclude yd = 0. 
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